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Arrays of coupled limit-cycle oscillators represent a paradigmatic example for studying synchro¬ 
nization and pattern formation. They are also of direct relevance in the context of currently emerg¬ 
ing experiments on nano- and optomechanical oscillator arrays. We find that the full dynamical 
equations for the phase dynamics of such an array go beyond previously studied Kuramoto-type 
equations. We analyze the evolution of the phase field in a two-dimensional array and obtain a 
“phase diagram” for the resulting stationary and non-stationary patterns. The possible observation 
in optomechanical arrays is discussed briefly. 
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Synchronization is an important concept in many 
branches of physics, chemistry, biology and other sci¬ 
ences [ ]. Within the past two years, a number of ex¬ 
periments have demonstrated for the first time synchro¬ 
nization between two nanomechanical oscillators [2-4]. 
These systems are driven through a Hopf bifurcation into 
a limit-cycle oscillation, where the energy pump is sup¬ 
plied through feedback or an optical drive. Future arrays 
of synchronized mechanical Hopf oscillators, with more 
than just two components, promise to provide robustness 
against both disorder and noise. Considerable theoreti¬ 
cal attention has recently been devoted to the problem of 
synchronization in arrays, both on the general level and 
for predicting the behavior of specific systems (e.g. in 
nanomechanics [5-10] or trapped ion systems [11]). Some 
progress has also been made in the quantum regime [12- 
16]. It is efficient to focus on the dynamics of the crucial 
phase degree of freedom, where the most prominent phe¬ 
nomenological model is the one introduced by Kuramoto 
[ .7, 18], which more recently has been supplemented by 
so-called reactive terms [6]. 

In the present work, we will explore synchroniza¬ 
tion and deterministic pattern formation [19] for a two- 
dimensional array of identical Hopf oscillators. We 
present the complete effective model for the phase dy¬ 
namics. Starting from the widely applicable model of 
coupled limit-cycle oscillators, we find that the classi¬ 
cal phase evolution is affected by extra contributions be¬ 
yond those investigated previously. These can have a 
significant impact on the dynamics. Our simulations of 
the effective model reveal various stationary and non¬ 
stationary patterns in different parameter regimes. Phase 
correlators, length scales, and macroscopic pattern dy¬ 
namics will be discussed. These are relevant for deter¬ 
mining whether an array can easily settle into a phase- 
locked state, which is important for applications. 

In the system that we study, all oscillators are un¬ 
dergoing motion on a limit cycle, see Fig. 1. We start 
with the following equations (“Hopf equations”), which 
describe the effective phase and amplitude dynamics of 



Figure 1: (a) A single Hopf oscillator undergoes dynamics on 
a limit cycle, with an amplitude set by external parameters 
(such as the drive providing the power) and a time-evolving 
phase (p(t ). (b) Array of coupled Hopf oscillators. Often, the 
system can be well described with a single degree of freedom 
per lattice site, the phase (fi . Due to the coupling, the phases 
can lock and phase patterns can form. 


these resonators close to the limit cycle [8, 13]: 

« = - 0 - (A, - i)g(i) + cos( Vl ), 

= + ( 1 ) 

Here O = Gt{A) is the frequency at the steady-state am¬ 
plitude A. The frequency-pull due to amplitude changes 
has been kept to leading order in the first line, 7 is the 
rate at which the amplitude is forced back to the limit 
cycle, and Fi(t) is the total force exerted on resonator i. 
For the physically most relevant spring-like coupling be¬ 
tween nearest neighbors, which will be studied here, the 
force is given as Fi = & X7j 7 Aj cos (pj). The coupling 
constant is k , and (j, i) indicates nearest neighbor sites. 

We consider the case of weak coupling k/(mQ 2 ) <C 1 
and assume 7/fi <C 1, (A/Cl)dfl/dA 1 . Then the 
amplitude fluctuations around the steady state value are 
small and the amplitude dynamics can be integrated out 
(for details on this step, also about disorder, see [8, 13] 
and the Supplemental Material (SM)). We arrive at effec¬ 
tive equations for the resonator phases. Keeping only the 
slow phase dynamics (and assuming identical resonators), 
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we get 

<Pi=C^2 cos (Vj - Vi) + Si ^2 sin (Vj - Vi) 

UA UA 

+ S 2 1 V V [sin( 2 ^ -<p k - <fi) ~ sin(^ fe - < p^] 
UA ( k A 

+ EE sm(ip k + ipj-2 (2) 

UA ( k A 

with C = k/2m£l, Si = (CA/j)(dSl/dA)\ A=A , S 2 = 
C 2 / 2 7 . We will call Eq. ( 2 ) the Hopf-Kuramoto model. It 
has been derived before in the context of optomechanics 
[ 8 , 13], but it holds generally for a set of weakly coupled 
Hopf oscillators. Our aim is to explore the dynamics of 
this model on a square lattice. 

The term sin(y?j — Pi) of Eq. (2) is well known from 
the Kuramoto model [17], or, equivalently, the XY 
model [20]. Here the term arises from the amplitude- 
dependence of the frequency [5]. Both contributions in 
the first line of Eq. (2) have been derived previously for 
coupled limit-cycle oscillators, see [21]. They are linear 
in the coupling k. In contrast, the prefactor S 2 is of sec¬ 
ond order in k. However, as discussed below, in realistic 
scenarios 7 and dVt/dA can become small, such that the 
regime S 2 ~ Si, (7 is easily reached. The S^-term can 
then have a profound influence on the pattern forma¬ 
tion dynamics. The additional contribution also displays 
next-to-nearest-neighbor coupling of the phases, in spite 
of the underlying intrinsic nearest-neighbor coupling in 
the lattice assumed here. 

We will first set the stage by highlighting several limit¬ 
ing cases of our model, some of which are known already. 
For (7 = 0, Eq. (2) can be rewritten in the form pi = 
—dU/dpi. Hence, the system will slide down to a mini¬ 
mum of the potential V . In contrast, for non-vanishing (7, 
such a potential does not exist and the system may never 
reach a stationary state (where pi is constant). The lim¬ 
iting case of Eq. ( 2 ) with S 2 = 0 has been studied before 
[22-26]. This is the Sakaguchi-Kuramoto model, usually 
written in the form pi—K ^ sin (pj — pi + a) with 
tan(a) = C/Si and K 2 = S 2 + (7 2 . 

The continuum limit of Eq. (2), which is valid for 
smooth phase fields, is given by 

p = SiAp - 2S 2 A 2 p - C(\7p) 2 + 4(7, (3) 

where Si = Si a 2 , £2 = S 2 a 4 , C = (7a 2 , with lattice 
constant a. In this model (with £2 = 0) it has been found 
that spirals can develop around singularities in the phase 
field [24, 27]. Besides, it has been analyzed in connection 
with chemical turbulence in one dimension [28-30]. 

The aim of this paper is to explore pattern formation in 
the full model, Eq. (2), in large two-dimensional arrays. 
Our main result is the pattern phase diagram discussed 
further below. The patterns we find will determine the 
phase synchronization dynamics of limit-cycle oscillators. 



Figure 2: Spiral patterns and length scales in the Hopf- 
Kuramoto model, (a) Stationary spiral pattern emerging from 
random initial conditions for Si/C = 2.0, S 2 /C = 0 on a 
N x N square lattice (N = 128) with periodic boundary con¬ 
ditions. (b) Vortex-anti-vortex pair (see inset) winding up to 
a stationary spiral-anti-spiral pair with a characteristic spiral 
arm width A. Parameters are like in (a), (c) Spatial cor¬ 

relations Re(exp(z(c^ m — p n ))) as a function of the distance 
\fm — r n | (rounded to the nearest integer). To obtain the 
data for (d), we extract the first correlation minimum posi¬ 
tion from parabolic fits and average over 10 runs with different 
random initial conditions, (d) The location of the first cor¬ 
relation minimum (red) and the spiral arm width A from (b) 
(black), as a function of the ratio Si/C, in units of the lattice 
constant. There can be hysteresis (blue). 


Our numerical results are mostly obtained from simu¬ 
lations with a random initial phase field, since that is the 
natural starting point in real systems (e.g. after switch¬ 
ing on the pump laser or the feedback driving the oscilla¬ 
tors into the limit cycle). After some transient dynamics, 
we often find patterns that do not change qualitatively 
any more on longer time scales. Moreover, in certain pa¬ 
rameter regimes, we find nontrivial stationary patterns. 

A typical final pattern of a simulation with large Si/C 
is shown in Fig. 2a. This pattern is stationary. It con¬ 
sists of many vortex-like “singularities”, where the phase 
changes by 2i r when going around in a closed loop. These 
points are surrounded by spiral structures. Spiral pat¬ 
terns in general are well-known as a recurring motif in 
pattern formation [19, 3f]. Since they form an important 
part of the patterns we observe, we now briefly discuss 
the properties of isolated spirals, produced from an initial 
condition with a vortex in the phase field (Fig. 2b). 

It is known that in related models, there is a tran¬ 
sition from stationary spirals to non-stationary spirals, 
i.e. a situation when the spiral centers are no longer 
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Figure 3: Pattern phase diagram of the Hopf-Kuramoto model, Eq. (2). Different colors indicate different patterns, which are 
discussed in the main text. Sharp transitions occur between stationary spirals and pulsating/mobile spirals (for small S 2 /C), 
and in the appearance of a 7r-defects”. Point markers indicate parameters explored by numerical simulation. Some typical phase 
patterns are shown in the insets (a) to (f). 


phase-locked to the bulk of the lattice [23, 32]. We have 
discovered that this transition also gives rise to a jump 
in the width of the spiral arms, A (Fig. 2). Outside of the 
jump, A increases with increasing S 2 /C and S\/C (black 
curve in Fig. 2d). When sweeping the parameter ratio 
Si/C up and down, we find hysteresis in the spiral arm 
width (blue line in Fig. 2d). The precise value at which 
the jump occurs can then depend on the parameter sweep 
rate. Our analysis illustrates that the microscopic details 
of the spiral center, on the scale of a few lattice sites, in¬ 
fluence both the spiral arm width and the macroscopic 
pattern considerably. Because the structure of the spi¬ 
ral core is complicated, we cannot provide an analytical 
prediction for A. 

We now turn to the statistical properties of the pat¬ 
terns which evolve out of random initial conditions (see 
Fig. 2a), as they are directly relevant for experiments. 
The spatial correlations of the phase field are character¬ 
ized by the correlator whose distance de¬ 

pendence is displayed in Fig. 2c. We find an oscillatory 
structure connected to the presence of spiral arms. On 
top of that, there is an exponential decay, due to the 
presence of many randomly located spiral centers. The 
position of the first minimum in the oscillations indicates 
the distance approximately set by half a spiral arm width. 
The dependence of this distance on the parameter S\/C 
is shown as the red line in Fig. 2d. Again, we find a sud¬ 
den jump associated with the transition from stationary 
to non-stationary spiral centers. We note that the spi¬ 


ral arm width A determined for isolated spirals does not 
agree completely with the length scale extracted from 
the oscillations of the correlator. The difference can be 
traced back to changes in the spiral core produced by the 
presence of other nearby spirals. 

There are only two dimensionless parameters, S\jC 
and S 2 /C , that determine the properties of the final 
pattern. Therefore, a complete overview of the various 
regimes in our model is provided by the “pattern phase 
diagram” in Fig. 3. This summarizes the main results of 
our studies, and we now explain its features. 

The transition discussed above, between stationary 
and non-stationary spirals, is sharp and can be traced 
up to intermediate values of S 2 /C. In addition, we find 
two classes of non-stationary spirals: “pulsating” spirals, 
where the core keeps orbiting in a small circle around a 
fixed location [32], and truly mobile spirals that move 
through the whole lattice. We will comment on their dy¬ 
namics later. We have not observed a sharp transition 
between the two regimes (Fig. 3). At larger S 2 /C, the 
transition is directly from stationary to mobile spirals. 

When decreasing the parameter Si/C even further, we 
find a crossover to “fluctuating” patterns, see Fig. 3c. 
These are non-stationary patterns with a complicated 
phase structure on the scale of the lattice. For the spe¬ 
cial case S 2 = 0, the location of the crossover (around 
Si/C ~ 0.8) matches the result obtained in [25]. 

The crucial macroscopic length scale of the Hopf- 
Kuramoto model, i.e. the spiral arm width, grows with 
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increasing S 2 /C. In view of that, it is surprising to see 
microscopic structures appearing at larger values of this 
parameter. Indeed, we find a sharp transition from the 
domain of “stationary spirals” to stationary patterns that 
contain “7r-defects”, see Fig. 3d. These are point defects 
which are offset by a phase difference of roughly 7 r from 
the smooth surrounding phase field. The stability of a 
single 7r-defect on a homogeneous background can be an¬ 
alyzed by semi-analytical linear stability analysis (in the 
limit (7 — y 0; for details, see the SM), which gives the crit¬ 
ical value S 2 /S 1 = 0.107. This defines the asymptote for 
the transition line in Fig. 3. Above the critical value, the 
7r-defect patterns form a fixed point of the dynamics and 
can be reached from random initial conditions. In con¬ 
trast to the pure spiral patterns, these patterns resolve 
the structure of the lattice and hence form a fundamen¬ 
tally different phase. Obviously, they cannot appear in 
the continuum model, Eq. (3). 

When increasing the parameter value S 2 /C further, 
the density of 7r-defects increases until we observe a 
smooth transition to “complex” patterns. These are sta¬ 
tionary patterns with a complicated phase structure on 
the scale of the lattice, see Fig. 3e. 

Finally, we note that the white region in the phase 
diagram could not be accessed due to the significant in¬ 
crease of timescales. Apart from that, we have discussed 
all phases in the Hopf-Kuramoto model, for positive pa¬ 
rameters. Changing the sign of C or Si will not give 
qualitatively different results: The emerging patterns can 
be reconstructed from the patterns discussed above by 
the transformations (frn,n —^ —< Pm,n for a sign change of 
(7, and (p m , n —tpm,n + (—l) m+n 7r/2 (“checkerboard 
gauge”) for a sign change of Si . This works for all values 
of S 2 . However, changing the sign of S 2 will lead to dif¬ 
ferent patterns. These involve structure on the scale of 
the lattice, where phase differences of roughly tt/2 play 
an important role. We will not discuss these patterns, 
because for coupled Hopf oscillators S 2 is positive. 

We now turn to a more detailed discussion of the spiral 
motion and interaction (see also [33] for the continuum 
case, at S 2 = 0). Whenever we observe mobile spirals, a 
fraction of the spirals and anti-spirals eventually annihi¬ 
late. In some cases, they can also be created dynamically. 
We observe that the spirals move through the array al¬ 
most independent of one another for small S 2 /C , whereas 
they tend to move in pairs for larger values of this pa¬ 
rameter. For large values of S 2 /C , mobile and stationary 
spirals can even coexist, see Fig. 4. Depending on initial 
conditions, the final state can then be non-stationary or 
stationary (if all mobile spirals annihilate). 

There is also a parameter regime where 7r-defects, sta¬ 
tionary and mobile spirals can all be present and inter¬ 
act: Upon the annihilation of a spiral-anti-spiral pair, a 
7r-defect can be left behind. This happens more often for 
larger values of S 2 /C. When a mobile spiral approaches 
a 7r-defect, it can induce the dissolution of the defect into 



Figure 4: Spiral motion in the Hopf-Kuramoto model, (a) 
Instantaneous phase velocity field (pi after a (long) inte¬ 
gration time T — 1000 /C. Mobile spiral centers are visi¬ 
ble as local inhomogeneities in Some stationary spirals 
(not visible) exist in the uniform regions. Parameters are 
S\/C = 1.6, S 2 /C = 0.5, N = 64. The corresponding phase 
field is shown in Fig. 3f. (b) Spiral positions at time T. The 
lines are the spiral trajectories (from T — 15/(7 to T; the tra¬ 
jectories have been slightly smoothened for clarity), (c) Spiral 
dynamics of a single spiral over a longer period of time. It 
can be seen that the spiral remains fixed for some time before 
starting to move again (this is usually induced via a kick by 
a nearby mobile spiral). The spiral preferably moves in the 
cartesian directions set by the lattice. 


a spiral-anti-spiral pair. However, the mobile spiral can 
also move across the defect and make it vanish. All these 
interactions play an important role even at late times. 

Experimental studies of the patterns discussed in this 
work could be implemented by direct local electrical read¬ 
out of the motion in electrically coupled nanomechanical 
resonator arrays [ ], or by optical readout of the motion 
in future optomechanical arrays [ 8 , 13] based on optome¬ 
chanical crystals [34, 35] or other platforms. The latter 
consist of an array of localized mechanical modes, each of 
them coupled to one localized optical mode, driven by a 
laser. The model parameters could be tuned by varying 
the laser power and detuning. Simulations of single op¬ 
tomechanical cells, where we extracted the phenomeno¬ 
logical parameters 7 , A and (dCl/dA)\ a=A> suggest that 
all the important regions of the pattern phase diagram 
could be explored. Near the Hopf bifurcation, 7 < C 
can be reached (since 7 0), so S2 > C. Furthermore, 

S 2 > Si holds as well for sufficient coupling fc, when 
A(dQ/dA) \ a= a < C. The motion can be read out by 
observing the light scattered from the sample. The in¬ 
tensity of the light scattered with wave vector transfer 
q is related to the structure factor (see SM), i.e. the 
spatial Fourier transform (at q) of the phase correlator 
7 t . In a typical experiment, the frequencies 
will be disordered, but first simulations (with standard 
deviation 0.1(7) do not show qualitative changes of the 
patterns we discussed. However, initially mobile spirals 
could be pinned at sites with lower frequencies [23]. 

The variety of patterns summarized in Fig. 3 are im¬ 
portant for synchronization dynamics and applications. 
For example, finite phase-differences across the array (in 











5 


stationary patterns) will reduce the total power output 
of the collective oscillator, and the mere presence of spi¬ 
rals can reduce the robustness against noise [9]. Finite 
frequency-differences (in non-stationary patterns) reduce 
the frequency stability. Tuning the parameters into suit¬ 
able regions will optimize the array’s properties. Future 
theoretical studies of the Hopf-Kuramoto model could 
include noise, which may lead to interesting effects, as 
discussed for similar models in [18]. In that context, as 
well as in the deterministic case, the role of spiral motion 
and interaction could be analyzed in more detail. 
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I. DERIVATION OF THE HOPF-KURAMOTO MODEL 


In this section, we derive the Hopf-Kuramoto model from the following general Hopf equations 


(pi — Eli 

Ai = -7 (Ai 


(A 

dA^= A[Ai 


Ai) + 


m 


— Ai) + 


m 

miEli(Ai) 


rfliEli (Ai)A' t 

sin (pi . 


COS (pi, 


(5.1) 

(5.2) 


Here, f \i{Af) is the amplitude-dependent frequency of the oscillator at site i, rrii is its mass, Ai is its steady-state 
amplitude and Eli = Eli(Ai) is the frequency at the steady-state amplitude. Other symbols have the same meaning 
as in the main text. The second term on the right-hand side of Eq. (S.l) arises from the expansion of Eli(Ai) around 
the steady-state amplitude Ai. For reasons that will become clear later, Eli(Ai) has not been expanded in the force 
terms in Eqs. (S.l) and (S.2). We assume that the oscillators are coupled by spring-like nearest-neighbor couplings 
so that the forces Fi(t) are given by 


Fi = kijXj = kijAj cos (pj, (S.3) 

UA UA 

where (j , i) denotes the nearest neighbors j of site i and kij = kji are spring constants. 

The derivation of the Hopf-Kuramoto model involves the adiabatic elimination of the amplitude fluctuations SAi = 
Ai — Ai , as well as leading-order expansions in the dimensionless, small parameters k^/ (ra^ 2 ), (. Ai/Eli)dEli/dAi 
and 7 /Eli. These parameters and the relative amplitude fluctuations SAi/A are assumed to be of the same order of 
smallness. Below, we will also only keep slowly varying terms. The derivation can also be found in [ ]. For some more 
details, see the Supplemental Material of [ ]. 

In order to eliminate the amplitude fluctuations, we rewrite Eq. (S. 2 ) in terms of the amplitude fluctuations and 
formally integrate the equation 

SA = - 7 SAi + % U* + 5A i) cos > (S.4) 

riflilLi yJ±i J ^ ^ 


to obtain the long-time limit result 

SAi(t) = — —7 a \ [ dt ' sinV % (Aj + 8Aj(t')) cos(S.5) 

miUi{Ai) J_ 00 

Since the integrand is proportional to to leading order it suffices to evaluate (pi(t f ) to zeroth order in the expansion 
parameters, i.e. (pi(t') ~ (pi(t) — Eli(t — t'). Thus, we find 

SAi ~ f ^ ^ An ((pi(t) — Eli(t — t')) cos — Elj(t — t ')). (S.6) 

(j i) m A^i\Ai) J -00 ' ' ' ' 

The integral can easily be evaluated. To leading order in ^/Eli, the result reduces to 

5A- ~ V Aikij Sin( ^ ~ ^ 

1 ~ miEli(Ai) 27 


(S.7) 
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To first order in the amplitude fluctuations, the equations of motion for the oscillator phases (S.l) can be expanded 
as 


Vi — £ 2 ^ 


dtti 

dAi 


| Ai=A^Ai + 


COS (fi 
mifii(Ai) 




COS ( Pj . 


(S.8) 


Corrections to ^i(Ai) ~ f^(A^) = are proportional to both (Ai/Q^dfli/dAi and 5Ai/Ai so that they are of second 
order. In the second term on the right-hand side, these are significant, but, since in the third term they are multiplied 
by another kij , they can be neglected there. Inserting (S.7), again replacing U^(A^) in the denominator by fb, we 
finally obtain 
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where we have also only kept slowly varying contributions by applying the approximations cos (pi cos cpj ~ E C os (cpi — 
( pj ), cos (pi cos pj sin (jpj —ph) — ^{sin(<^ ~Vk) — sin (vi + Vk — 2 Vj)} and cos <pi cos <pj sin(^ — Lp k ) — |{sin(2 Vi ~Vj ~ 
2 (p k ) — sin (<p k — Vj)}' I n the special case of identical oscillators = H, rrii = m and Ai — A for all i and uniform 
couplings k^ = k for all neighbors i, j, this obviously reduces to equation (2) in the main text, where we have also 
neglected the trivial term — Cl on the right hand side. 


II. SEMI-ANALYTICAL STABILITY ANALYSIS OF POINT DEFECTS 

In this section, we present the analysis of the stability of a single 7r-defect on a homogeneous background phase 
field for the Hopf-Kuramoto model (Eq. (2) in the main text) with (7 = 0. For this case, the aforementioned phase 
configuration, which we call , is a fixed point of the dynamics, i.e. <p 9 = 0 for all sites i. Besides, the equation of 
motion can be written as ipi = with the potential 

U(if(p N 2) = El Ev(l - cos(^j - <p*)) + s 2 [ sin(<£,- -Pi)] }■ 

i CM) CM) 

We calculate the Hessian d 2 U/{dpi&Vj) and evaluate its eigenvalues for the phase configuration cp° numerically. If 
at least one of the eigenvalues is negative, the configuration is unstable. A single eigenvalue, corresponding to the 
translational mode of the system, might vanish without disturbing our analysis. We always find this zero eigenvalue. 
For small values of S 1 /S 2 , all the other eigenvalues are positive, which means that 7r-defects are stable (see Fig. 
S.l). With increasing S 1 /S 2 , the eigenvalues change linearly with this parameter. A single eigenvalue A _ (Si/S 2 ) has 
a negative slope, so it becomes negative at some critical value (S 1 /S 2 Y ~ 9.34, rendering the phase configuration 
unstable. This gives the (inverse) value S 2 /S 1 ~ 0.107 given in the main text. 


III. READ-OUT OF THE MECHANICAL RESONATOR PHASE FIELD 

Here, we will show how the intensity of the light reflected from an optomechanical array is related to the spatial 
Fourier transform of the phase correlator. The intensity of the light reflected from an optomechanical array with 
lattice sites at Tj is given as 

| J E , (^)| 2 /|.E ,in | 2 = ye-«e^| 2 . (S.10) 

3 

The phase of the light reflected from site j is 0j = # max cos(y?j), where (9 max depends on the system parameters. 
If the mechanical frequency is much smaller than the cavity intensity decay rate ft, then # max = AG/ft, with the 
mechanical amplitude A and the optical frequency shift per displacement G [3]. For small # max , Eq. (S.10) can be 
expanded and we get 
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Figure S.l: The eigenvalues A m of the Hessian d 2 U /(dtpidipj) for a homogeneous phase configuration with a single 7r-defect, as 
a function of the parameter S 1 /S 2 for a lattice of size 20x20. Note the eigenvalue A - leading to the instability of the 7r-defect. 
The zero eigenvalue X T belongs to the translational mode of the phase configuration. 



Figure S.2: (a) Partial intensity e i ^' rj 0j\ 2 j of the light reflected from the phase field given in Fig. 2a in the main text. 

Parameter 0 max — 0.01. (b) Partial intensity e l ^' rj 6j | 2 ^ in dependance on the radial coordinate q r — y/q% + , averaged 

over 11 different random initial phase configurations (black) and for a single phase configuration as in (a) (blue). 


We average over time, use (9j)t = 0 and (9j)t = (# max ) 2 /2 , and arrive at 

/ \ (f) max\2 . . 

(\E(q)\ 2 /\E™\ 2 ) t = (1 - V_^)|^ e ^| 2 + (Ee^9j| J ). (S.12) 

3 3 

For large arrays, the first term will only give contributions very close to q = 0. For small arrays, these contributions 
may have to be eliminated by calibrating the measurement device with a known phase field. The second term can be 
evaluated to give 


, , /nmax\2 

(|^y<?^0.| 2 ) = l _ L ^ e -*9-(n-^) Re ( e i(w-^)) t . (S. 13 ) 


On the right-hand side of this equation, the discrete Fourier transform of the correlations in the system appears. 
We have analyzed similar correlation functions in connection with the spiral length scale, see Fig. 2 in the main text. 















9 


From Eq. (S.13) we see that we can learn about the correlations by detecting the intensity of the reflected light. An 
example for the part of the detected light intensity that is given in Eq. (S.13) is given in Fig. S.2. 
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